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Abstract 

Xcj (J=0,l,2) production associated with a charm quark pair in hadron collisions is calculated in 
the framework of nonrelativistic QCD at the Tevatron and LHC. It is found that the color-singlet 
contribution is small and the color-octet contribution may be dominant in the large px region. 
The differential cross section of pp — > Xc + cc is at least one order of magnitude smaller than the 
next-to-leading order result of XcJ inclusive production pp — > Xc + X, therefore XcJ production 
in pp — > Xc + cc may have negligible influence on the ratio R Xc = measured by CDF at the 
Tevatron. The feeddown contribution from Xcj + c + c to J /tp + c + c is found to be large compared 
with J ftp direct production and may have important influence on the measurement of J ftp + c + c. 
The validity of fragmentation approximation for the process is also discussed. 

PACS numbers: 12.38.Bx, 13.60.Le, 14.40.Pq 
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I. INTRODUCTION 



Charmonium production associated with a cc pair is a good experimental observable in 
understanding the production mechanism of heavy quarkonium. The associated production 
has been extensive studied in the literature. In e + e~ annihilation at B factories, J/ip + c + c 
was found to have a very large fraction of J/ip inclusive production Jl], 0. This phenomena 
can be understood by a large next-to-leading order (NLO) QCD correction to the color- 
singlet (CS) J/ift + c + c pr o duct ion 0, and a relatively small NLO QCD correction to 
the CS J/ip + X non _ C c pro duct ion ||. These studies also imply that the color-octet (CO) 
contribution to J/ip production may be very small and even negligible, and a set of severe 
constraint on the linear combination of related CO matrix elements was further obtained by 
analyzing J/ijj production in e + e~ annihilation ||. At LEP, in the Z° decay, J/^+c+cp was 
found to be the dominant contribution to the J/ip inclusive production^. On the contrary, 
in 77 collisions, the contribution of J / ip associated production || was estimated to be several 
orders of magnitude smaller than the experimental data of J/ip inclusive pro duct ion ||. In 
hadron collisions at the Tevatron, theoretical predictions |T0|, [TT|j showed that the J/ip + c + c 
contribution was significant in the large pt region compared with the NLO result of non- 
ce contributions, and the produced J/ip is mainly unpolarized, which is analogous to the 
polarization of J/ip in inclusive production ]I2], [HJ. The integrated cross section was also 
significant and showed a great measurable potential both at the Tevatron and RHIC. J/ip 
associated production was also considered in the T decay [14[] to explore the CO mechanism 
in heavy quarkonium decays. As shown, most of those studies focused on J/-0 associated 
production. However, due to the importance of charmonium associated production, studies 
of associated production for other charmonium states other than J/ip may also be valuable. 

In this paper, we perform a calculation for the P-wave charmonium Xcj (J=0,l,2) associ- 
ated production in hadron collisions in the framework of nonrelativistic QCD(NRQCD) fllS] . 



The motivation for this work is two fold. The first is related to the ratio i? v = measured 
by the CDF collaboration at the Tevatron. CDF found that R Xc approaches to about 0.75 
at large pr|[6|]- However, if the Xcj inclusive production is dominated by the CO process 
(because the CS contribution at large pt is suppressed by ■jf at leading order(LO)), the 
value of R Xc should tend to be |, which is predicted by naive spin counting. Recently, the 



calculation of NLO QCD correction to Xcj + X non _ c5 inclusive production is performed |T7 
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and it is found that the NLO correction for CS channel can bring out a -4- term, which 
makes the CS contribution much important at large pt, and then by combining CS with 
CO contributions one is able to fit the experimental value of R Xc quite well over a wide pt 
region. The Xcj associated production with a cc pair is of the same order in pQCD as the 
NLO Xcj + X inclusive production and it contains also fragmentation contributions which 
scale as So it is interesting to check whether the Xcjcc associated production is also very 
large and whether it can further improve theoretical predictions of ratio R Xc . The second 
reason is that the measurements at the Tevatron for the production rates of J/ip + c + c 
and J/i/j + X are important on shedding light on understanding the J / ifj production mech- 
anism in hadron collisions. And the prompt J/ip production receives significant feeddown 
contributions from ip(2S) and Xcj- So it is important to know how large is the feeddown 
contribution to J/i/j associated production from Xcj associated production. The result itself 
in this work can also give information for directly measuring Xcj from Xcj + c + c production 
at hadron colliders. 

The remainder of this paper is organized as follows. In Sec. |I| we briefly describe our 
calculation method. In Sec. | we give our numerical results and analyze the obtained 
results. In the last section, we give a summary. 



II. CALCULATION OF pp — >Xcj + c + c 

In the framework of NRQCD factorization, the cross section for the Xcjcc associated 
production in proton-antiproton collisions has the following form: 



da[pp Xcj + c + c] = ]P J dx 1 dx 2 fi/ p (x 1 )f j/ p(x 2 ) 

i,j,n 



x da(i + j -> cc[n] + c + c)(0 XcJ [n}) (1) 

where i, j denote the initial state partons from the proton or anti-proton. We assume the 
contribution from light quark annihilation to be negligible, so i, j are gluons in our case. The 
quantum numbers n represent the color and orbital angular momentum of the intermediate 
cc states at short-distances, which evolve into the Xcj meson at long-distances. At leading 
order in relative velocity v of the cc pair, n can be taken as n = 3 P [ j ] for the CS and n = 3 Sf 
for the CO intermediate states. We use FeynArts[lS] to generate Feynman diagrams. For 
the CS case, there are 56 Feynman diagrams, of which some representative are shown in Fig. 



|]. For the CO case, there are 16 extra Feynman diagrams relative to the CS case and their 
extra topology structures are shown in Fig. 0. These extra topology structures actually 
represent gluon fragmentation contributions. 




(m) (n) 



FIG. 1: Representative Feynman diagrams for pp — > Xcj + c + c. 




(e) 



FIG. 2: Representative extra diagrams for CO channel. 

The fragmentation diagrams can give a contribution for the differential cross section 
d J^s , while remained diagrams can at most give ^- terms. Therefore, one may expect the 
fragmentation contribution is dominant at large p?- Note that, however, not only the CO 
channel ( 3 iS{ 8 ^) have fragmentation contribution, but the CS channel has also quark (anti- 
quark) fragmentation contributions. In Fig.|l], the diagrams (b) and (k) represent (anti- 
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) quark fragmentation contributions in the large px approximation, where one of the charm 
or anti-charm quark fragments to Xc plus another charm or anti-charm quark. However, only 
when pt is large enough, these quark diagrams can reach the fragmentation region and give 
leading contributions. Otherwise, these quark diagrams can give only small contributions. 
We will discuss the fragmentation approximation in the next section. 

We use the spinor-helicity method to deal with Feynman amplitudes ||19|| and use the 



package S@M||20| to simplify the Feynman amplitudes in spinor form. Based on this method, 



the spinor form for massive external fermions can be written as 

«±§(P) = ^^==0* + m)|g ±>, (2) 

where u ± i(p)(v±i(p)) are Dirac spinors of massive (anti-)fermion with momentum p and 
spin ±|; qo± are reference Weyl spinors with light-like reference momentum go and helicities 
A = ±1. In principle, "cfo should be along the axis of to guarantee the validity of 
above equations for individual spin. However, if we just concern about a result by summing 
over spin, go can be chosen arbitrarily. Here, it is chosen to be one of the initial partons' 
momentum in order to simplify the calculation. The polarization vectors for external gluon 
with momentum k and light-like reference momentum go are represented as 

t + (k,q ) = 7 ^r(|fc-)<?o-| + Ifc+XM), ( 4 ) 

{qo\k} 

f(k,q ) =j£^(\k + )(q 0+ \ + \q -)(k.\). (5) 



Using the following identity, we decouple the spin projection operator |2T[ for bound states 
from Feynman amplitudes: 

~ m c)\lox 2 }{<lox2\ p isJqo\ s }{qo\ s \(f c ~ m c ) = 2p c ■ q 2p 5 ■ q Q P 1Sz , (6) 

A2A3 

where there is a relative velocity v between p c and p 5 for the P-wave case. Then with the 
help of Fierz transformation (and its generalized forms), the amplitudes are reduced to 

Mi = CfViA, (7) 
where i is the index for different diagrams and /,■ are the simplified fermion chains, and the 



three-gluon vertex are properly dealt with (see [19] for details). The specific expressions of 
fj encountered here are listed in the appendix. 



We write the polarization tensors for Xcj explicitly. For Xco it is symmetric for the two 
indexes \w: 

<T = + ^; (8) 
for Xci it is anti-symmetric for the two indexes fiu: 

(9) 



V2' 
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for Xc2 it is again symmetric for the two indexes \xv: 
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The definition of e x , e y , and e z are 



(10) 



sf=^(e^-e^); (11) 



(12) 
(13) 
(14) 
(15) 
(16) 



e x (P) = (0, cos 9 cos 0, cos 6 sin 0, — sin#), (17) 
e y {P) =(0, - sin 0, cos 0,0), (18) 
e z (P) =— (|^|,Posm0cos0,Posin0sin0,PoCOse), (19) 



where M, p\ and P u are the mass, momentum and energy of Xcj', angles 6 and describe 



Xcj's direct ion p2|. For CO 3 Sf\ the spinor-helicity forms of polarization vectors are kept 
as (<7oa|^|(Zoa)) (<?oa|^|<?oa), (<Zoa|W|(Zoa) until numerically squaring the amplitudes in Fortran 
program. 

In our numerical calculation, we give only a rapidity cut condition for \ c . However, 
to detect the associated production, one should detect at least another hadron containing 
charm or anti-charm quark, and the rapidity cuts from experimental facility should also 
apply for the (anti-)charm quarks in principle. 
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For phase space integration, we use the general 2 — > 3 phase space expression, plus two 
fold momentum fraction integration for initial partons: 

-T- = / dzi/^Cxx) / dx 2 / 3/ ,(x 2 ) / dfc 5 ° / dk° 3 / dr/V|M| 2 , (20) 

a PT J5 JS/xi JyJml+ P 2 T J(k°) mi n Jo 

where 5 = 16m + Ap r ^ j g ener gy Q f ^ j g ener gy Q f one Q f th e emitted charm or 
anti-charm quark, rj describes the angle between the plane for the final three particles and 
the plane chosen for the two initial partons, and we omit flux factor and other normalization 
factors. The upper and low limits for k® and k® integration are a little complicated so 



we don't list them here. We use Vegas |23[ in Fortran program to perform the numerical 
integration. The correctness of our phase space integration program is verified by comparing 
the calculated J/ip + c + c production (we calculate it again) with the result from Ref . 1 10| . 



III. RESULT AND ANALYSIS 

In numerical calculation, we choose m c = 1.5GeV. The factorization scale and renormal- 
ization scale are both chosen as /j,q = ^p 2 T + Am 2 c . We use CTEQ6M as PDF input. The CS 
matrix element (0 XcJ [ 3 Pj 1 ']) is related to the P-wave function at the origin by the formula: 
(0^ J [ 3 P [ j ] ]) = {2J+ £ 3Nc \R' P (0)\ 2 and we choose: \R' P (0)\ 2 = 0.075GeV 5 from the potential 
model calculations |24]. For the CO matrix element (0 XcJ [ 3 S'{ 8 ']), by spin symmetry, we have 
the following relation that (0*<*[ 3 Sf ] ]) : (0* cl [ 3 Sf ] ]) : (0* c2 [ 3 Sf ] ]} = 1 : 3 : 5 and we use 
(O Xc0 [ 3 S'{ 8 ']) Rj 2.2 x 10~ 3 GeV 3 as the central value obtained from fitting the Xcj inclusive 



production pp to Xcj + X non _ cS at the Tevatron|T7j|. 

In Fig. |^, we show both CS and CO contributions to the differential cross section of 
Xc + c + c. It is found that the CO contribution dominates over production in the large 
Pt region, and it decreases much slower than that of CS as pt increases. This seems to 
conflict with that both of the two channels should behave as ^- at large px- To see this 
point more clearly, we fit the parton differential cross section (taking out the influence of 
PDF) with and it turns out that the CS channel scales roughly as while the CO 
channel scales as -\- in the region px < 17 GeV. This implies that the CO contribution in 
Fig. ^ is dominated by gluon fragmentation just as expected, but the CS channel has not 
reached the (anti-) quark fragmentation region, and its contribution is still suppressed in the 
moderately large region (e.g., px < 17 GeV). 
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FIG. 3: Differential cross sections of XcJ+c+c as functions of transverse momentum at the Tevatron 
with y/s = 1.96GeV and rapidity cut \y XcJ \ < 0.6. The dashed line denotes CO contribution, and 
the solid line denotes CS contribution. 



S 



In order to understand the above mentioned px behavior, we study the J/ip+cc production 
as an example. We calculate (anti-)quark fragmentation diagrams in axial gauge to include 
all 4p contributions, yfs can be set to be lOOTeV to enable us to calculate at px as large 

Pt 

as possible. We find that when px = 50GeV, (anti-)quark fragmentation contribution has 
a fraction of about 70% of the total differential cross section, and then the fraction rises 
to about 93% when px = 150GeV. And the fraction reaches 100% (within the calculation 
errors) when px > 450GeV. Thus we find that (anti-) quark fragmentation approximation 
is only valid for very large pt {px > lOOGeV at least), while for the presently interested px 
region (px < 17GeV) the \- term induced by (anti-)quark fragmentation is very small and 
not important. This explains the fact that the CS channels in Xcj production behave almost 
as -i . 

In Fig. [|, we depict the differential cross section for Xc + c + c and the NLO result for 
Xc + X non _ C c [0 as comparison. We find that the contribution from Xc + c + c is about 



two orders of magnitude smaller than Xc + A non _ cE at small px- The fraction of Xc + c + c 
in total Xc + A increases gradually and reaches at most 20% at px as large as 60GeV. The 
smallness of the pp — > Xc + c + c fraction lies in the fact that one of the main sources of 
Xc + X non _ cc - is the CO, which scales as 4r and begins its contribution at order a^, while 
the dominant contribution of Xc + c + c is suppressed by both a s and phase space. Based 
on this analysis, we may conclude that Xc + c + c has negligible influence on the Xc inclusive 
production. As a result, fitting the ratio R Xc measured by CDF|FB| can not be improved by 



including pp — > Xc + c + c as compared to the pp — > Xc + A raon _ cg result [17j • 

In the NLO prediction [|17|, the feeddown contribution of Xcj + A non _ cg possesses about 
30% of the prompt J/tp production rates at px = 20GeV at the Tevatron, and it can 
give a great influence on J /if; prompt production. Thus we also evaluate the feeddown 
contribution of Xcj to J/ip + c + c to see whether this contribution is also large. In the 
calculation, we ignore the difference between px of J/if> and Xc- Note that the feeddown 
from Xcj may have important influence on prompt J/^'s polarization. This effect relies on 
Xcj's polarized production rates and also the helicity amplitudes of Xc radiative decays to 
polarized J/if>. The related formula can be found in Ref.[^]. In this work, we only consider 
the unpolarized Xc + cc production but ignore the polarization effects. The branching ratios 
for Xcj radiative decays to J/if> are Br(xco — > J/i-' + l) ~ 0.013, Br(xd J/i' + l) ~ 0.36, 
Br(xc2 ~> J/i> + 7) ~ 0.20 respectively |2B[]. In Fig. ^, we give the feeddown contribution 
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FIG. 4: Comparison of the differential cross sections of XcJ + c + c with the NLO results of 
Xcj + ^non~cc at the Tevatron with y/s = 1.96GeV and rapidity cut \y Xc \ < !• 
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FIG. 5: Contribution of Xcj feeddown to prompt J/ip + c + c + X production at the Tevatron with 
rapidity cut |y Xc | < 0.6. The dotted line denotes the contribution from CS XcJ feeddown, the short 
dashed line denotes that from CO Xcj feeddown, the long dashed line denotes that from CO+CS 
XcJ feeddown, and the solid line is the direct J/ip + c + c contribution (from (|lO|l). 

from XcJ as a function of px- By comparing it with the direct production, we can see that 
the feeddown contribution from Xcj is small in the low p? region, but it is about a factor 
of 2 greater than direct J/ip + c + c contribution when > 15GeV. The turning point is 
at pt ~ 9GeV where the feeddown contribution begins to exceed the direct contribution. 
For Xc + c + c is dominated by the 3 S'{ 8 ' channel at large pr, one may anticipate that 
the CO J/ip + c + c contribution from 3 S'{ 8 ' may also play an important role in the direct 
J/ip + c + c production. However, the magnitude depends on the size of (0 J ^( 3 Sf^)}. In 
a recent work ||27|| , the authors find (0 J ^( 3 Sf'^)) might be small. As a result, Xc feeddown 
could become the main source for prompt production of J/ip + c + c, if (0 J ^( 3 S^)) is 
small. Thus when measuring the production cross sections for prompt J/tp + c + c at hadron 
colliders, the feeddown effect from Xc + c + c can be very important and should be taken into 
consideration. We also note that this situation is different from that at B factories, where the 
Xcj + c+c production rates in e + e~ annihilation for both CS and CO are very small, and their 
feeddown contributions to J/ip are negligible, therefore, the NLO J/ip + c + c theoretical 
results (including direct and ip(2S) feeddown contributions) are basically consistent with 
experimental production rates. At hadron colliders, however, the Xcj + c + c feeddown effect 
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becomes more important. 

We also give the prediction of Xcj + c + c associated production at the LHC with y/s = 
7TeV. For the CMS detector, the rapidity cut is \y Xc j\ < 2.4 and for the LHCb detector, the 
rapidity cut is chosen as 2 < y Xc j < 4.5. The results are shown in Fig. |6] and Fig. ^ including 
the CO contribution, CS contribution and the total differential cross sections. We find that 
the differential cross sections for Xcj + c + c associated production at the LHC show a similar 
behavior as that at the Tevatron: the CO contribution is much larger than CS contribution 
in almost all pr region. So Xcj associated production cross section is dominated by the CO 
contribution. As a result, measuring Xcj + c + c production can be used to determine the 
CO matrix element (O Xc0 [ 3 5'{^]}. Predictions for LHC with \/s = 14TeV and \y Xc j\ < 3 are 
shown in Fig. |8|. 



Summary 

In this paper, we investigate the Xcj associated production with a charm quark pair 
PP ~~ >* Xcj + c + c at hadron colliders in the framework of NRQCD. By comparing the 
differential cross sections i.e. the transverse momentum distributions for Xcj associated 
production, we find that CO dominates the production rate at large px- Also, the differential 
cross sections for associated Xcj production are at least one order of magnitude smaller than 
the NLO result for Xcj + X non _ c5 . As a result, Xcj + c + c production has negligible influence 
on the R Xc value measured by the CDF collaboration. We also evaluate the Xcj feeddown 
contribution to prompt J ftp + c + c production, and find that the feeddown contribution 



is very large compared to direct J/ip + c + c pro duct ion |L0| at large pr, which illustrates 
the importance of Xc feeddown effect in the measurement for J ftp associated production 
cross sections and polarization parameters. The fragmentation approximation is analyzed 
and our conclusion is that the fragmentation contribution is dominant for the CO channel, 
while for the CS channel the fragmentation diagrams' contribution dominates over the total 
differential cross section only at the px > lOOGeV region. 

Finally, we note that in the Xcj associated production, the LO result in a s has already 
contained the 1/pr 4 term, which is the leading term in 1/pt expansion at large pt, and high 
order corrections in a s can at most give the 1/pr 4 term, but suffer from suppressions due 
to extra powers of a s . So we expect that high order corrections in a s can not significantly 
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FIG. 6: Differential cross sections of Xcj + c + c as functions of transverse momentum at the LHC 
with y/s = 7TeV and rapidity cut |y XcJ | < 2.4. The dashed line denotes CO contribution, and the 
solid line denotes CS contribution. 
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FIG. 7: Differential cross sections of Xcj + c + c as functions of transverse momentum at the LHC 
with \fs = 7TeV and rapidity cut 2 < y XcJ < 4.5. The dashed line denotes CO contribution, and 
the solid line denotes CS contribution. 
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FIG. 8: Differential cross sections of Xcj + c + c as functions of transverse momentum at the LHC 
with y/s = 14TeV and rapidity cut |y XcJ | < 3.0. The dashed line denotes CO contribution, and 
the solid line denotes CS contribution. 
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change the results obtained in this work. Another notable result in this work is that the CO 
contribution dominates over Xcj + c + c production at large pr (say, px > 7GeV), therefore 
measuring the process at hadron colliders, especially at the LHC, may provide important 
information for the production mechanism of heavy quarkonium, while the experiment itself 
may also be very interesting and challenging in view of the complexity of the measurement. 
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Appendix 

We list all the fermion chains encountered in our calculation as follows. 

/o(gi,?2, Ai, A 2 ) =(goA 1 |(^ 1 + rn)y 5 ($ 2 -m)\q 0X2 ), (21) 

<?2, k, Ai, A 2 , A 3 ) =(goAil(j/ 1 + rn)y s (fi -$ 2 + m)£ X3 (k, g )(^ 2 - rn)\q x 2 ), (22) 

f2(qi,Q2,k, Ai, A 2 , A 3 ) =(g 0Al |(^ +m)t X[i (k,q )($ 1 - # + m)-f S (<j 2 - m)\q 0X2 ), (23) 

f 3 (q u ? 2 , k, A l5 A 2 , A 3 ) HcoaiK^ + ™)i M {i 2 - ™)koA 2 ), (24) 
f4(Qi,Q2,kx, k 2 , Ai, A 2 , A 3 , A 4 ) =(g 0Al |(^ + m)y s (h + h - $ 2 + m){ X3 (k 1 ,q ) 

(h ~$ 2 + rn)^ 4 (k 2 , q )(^ 2 - m)\q 0X2 ), (25) 
fs(qi,q2,fa, k 2 , Ai, A 2 , A 3 , A 4 ) = (qo\ 1 \($ 1 + m)^ 3 (k 1 ,q )(^ 1 - & + m)y s 

(h ~ $ 2 + m)^ X4 (k 2 , q )(q- 2 - m)\q 0X2 ), (26) 
h{qiA2,fa,k 2 , Ai, A 2 , A 3 , A 4 ) =<g Ai|(^ 1 + m)^ X3 (k u q )($ 1 - & + m)^ 4 (/c 2 , q ) 

(^ - ^ - % 2 -m)7 5 (^ 2 -m)\q 0X2 ), (27) 

frfai,®, fa, k 2 , Ai, A 2 , A 3 , A 4 ) =(qo\ 1 \($ 1 +m)^ Xs (k 1 ,q )^ Xi (k 2 ,q )'ys(^ 2 -'m)\qox 2 ), (28) 

h(qi,q2,fa,k 2 , Ai, A 2 , A 3 , A 4 ) =(g A 1 |(j# 1 + m^s^ifa, q )^ X4 (k 2 , q )($ 2 - m)|g A 2 ), (29) 

/g(?i,?2, fci,fc2, Ai, A 2 , A 3 , A 4 ) =(g 0Al |(^ +m)/ A3 (/;; 1 ,g )(^ 2 -m)|g A 2 ), (30) 
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where qf = q\ = ml, k ,ki,kz and qo are light-like vectors. 
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